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Many non-equilibrium processes on scale-free networks present anomalous critical behavior that is 
not explained by standard mean-field theories. We propose a systematic method to derive stochastic 
equations for mean-field order parameters that implicitly account for the degree heterogeneity. The 
method is used to correctly predict the dynamical critical behavior of some binary spin models 
and reaction-diffusion processes. The validity of our non-equilibrium theory is furtherly supported 
by showing its relation with the generalized Landau theory of equilibrium critical phenomena on 
networks. 



The discovery that many real systems can be repre- 
sented, at some large-scale level, as complex networks 
with a disordered topology and an highly heterogeneous 
structure, has recently motivated a genuine interest for 
studying the properties of dynamical processes occurring 
on top of them [l| . The primary statistical feature that 
seems to have a strong influence on dynamical proper- 
ties is the existence of a very heterogeneous spectrum 
of local interaction patterns, reflected at the topological 
level into a very broad degree distribution. A deeper un- 
derstanding of its influence on the dynamics is obtained 
isolating this feature from the remaining complexity of 
real networks by considering ensembles of heterogeneous 
random graphs with power-law degree distribution . 

Despite the local tree-like topology, random graphs 
present long loops ensuring that the average distance be- 
tween any two points of the system scales at most loga- 
rithmically with the system's size Therefore, physical 
properties on random graphs are not explicitly affected 
by spatial fluctuations. Naively, this suggests that, in the 
absence of additional disorder and frustration, the gen- 
eral behavior of equilibrium and non-equilibrium prop- 
erties on the ensemble of uncorrelated random graphs 
should be of mean-field (MF) type. Nevertheless, analyt- 
ical and numerical studies of statistical mechanics mod- 
els defined on scale-free networks have revealed that fat- 
tails in the degree distribution can generate non-trivial 
critical behaviors in correspondence of continuous phase 
transitions 0]. In networks with degree distribution 
P{k) ~ k-^, all moments Zp = (kP) = Y^k^^Pi.^) with 
p > 7 — 1 diverge with the system's size N; therefore 
scale-free networks with 2 < 7 < 3 present a finite av- 
erage connectivity but uncontrolled degree fluctuations. 
On infinite networks, diverging degree fluctuations can 
prevent the occurrence of phase transitions, and when- 
ever they exist the position of critical points and critical 
exponents explicitly depend on the value of 7. 

The standard technique to investigate these phenom- 
ena on random graphs is the heterogeneous mean-field 
theory (HMF) 4], namely by deriving a system of cou- 
pled mean-field equations for degree-dependent quanti- 
ties. The approach is based on the annealed network 



approximation, in which network's topological character- 
ization is strictly statistical, given in terms of the degree 
distribution P{k) and, possibly, multipoint degree cor- 
relations. Once the system of equations is solved, one 
computes the average global behavior in a self-consistent 
way to obtain the 7-dependent anomalous MF critical 
behavior. A complementary approach proposed by Golt- 
sev et al. [1| consists in assuming some reasonable form 
for the effective free-energy functional of the system and 
deriving from it the critical behavior of the corresponding 
MF order parameter. This method generalizes the mean- 
field Landau- Ginzburg theory to take into account the 
topological heterogeneity of the system and provides a 
general conceptual framework to understand equilibrium 
critical phenomena on networks. A similar general the- 
ory for non-equilibrium processes around critical points is 
still lacking, even if an explicit dependence of the critical 
scaling behavior on the degree distribution's exponent 
7 has been observed in several non-equilibrium models 
by means of the HMF theory and numerical simulations 
[3 0, B B U, ^ main progress was recently 

achieved by showing that, like in equilibrium, the finite- 
size scaling does not depend only on the system's size TV, 
but it may also depend on the upper cut-off kc of the 
degree distribution [l^, [13] ■ The non-equilibrium criti- 
cal dynamics is therefore crucially influenced by the tail 
of the degree distribution and by the behavior of kc{N), 
particularly in the region 2 < 7 < 3. 

In this Letter we provide a general mean-field ap- 
proach taking into account the anomalous effects intro- 
duced by degree heterogeneities, that could be applied to 
many classes of non-equilibrium processes in a simple and 
systematic way, generalizing the equilibrium Landau- 
Ginzburg theory proposed by Goltsev et al. We will 
first expose the general theory, for uncorrelated random 
graphs in the annealed network approximation, then its 
validity is verified by direct application to three proto- 
typical models of non-equilibrium dynamics on networks: 
the voter model, the Ising-Glauber dynamics and the 
reaction-diffusion processes. We show that degree hetero- 
geneity affects the non-equilibrium dynamics in different 
ways depending whether or not the system's size N and 
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the cut-off kc are finite. In networks with a finite maxi- 
mum degree, we eventually recover the mean-field scaling 
laws by tuning the control parameters sufficiently close 
to criticality. On the contrary, in networks with a cut-off 
rapidly diverging with the system size, the dynamics in 
the thermodynamic limit is governed by universal scaling 
laws that depend explicitly on the exponent 7. 

General Theory - In a general dynamical process on 
network, the microscopic dynamics is described by vari- 
ables {xi} defined on the nodes i — 1, . . . ,N; we con- 
sider their degree-dependent mean-field averages Xk ~ 
Z^i:fe,=fc Si '^fei.fc- The corresponding HMF rate equa- 
tions are of the form Xf^ = /(X^, X), where X is a global 
average quantity to be determined in a self-consistent 
way averaging over the degree distribution. 
In some cases the dynamics of X is slow compared to the 
microscopic time-scales. This is evident i) when X is a 
conserved quantity (e.g. in the voter model), but it is ex- 
pected to be correct also ii) when the microscopic dynam- 
ics is driven by the behavior of high-degree nodes (e.g. for 
reaction-diffusion systems and Glauber dynamics) [2l| . 
If X is slow with respect to the relevant microscopic dy- 
namics, we can use a quasi-static approximation impos- 
ing that Xk rapidly relaxes to a quasi-stationary state 
that depends parametrically on the global quantity X{t). 
The condition X^ — allows to express Xk = g{k, X{t)) 
as a function of k and X{t). 

The next step consists in deriving a stochastic equation 
for the slow variable X{t). The mean E[AX] and the 
variance V[AX] of the variation of X during a unit time 
interval can be easily computed. More precisely, one has 
to consider the variations of X that are due to single 
microscopic update events on nodes with degree fc, then 
replace the expression Xk = g{k, X{t)) and average over 
the degree distribution. At long times, close to the crit- 
ical point, X is small, so we can expand these terms in 
powers of X obtaining a Langevin equation of the type 
X{t) ^ aX + bX^ + ■■■ + y^{a'X + b'X^ + ...)/Nr){t), 
where r]{t) is a gaussian noise. Because of the expansion 
of g{k, X) for small X, all coefficients are functions of the 
moments Zp of the degree distribution. In the N 00 
limit, noise terms will drop out but some of the other 
coefficients a,b, . . . might diverge. 

One can use an action formalism to verify that a Langevin 
equation with diverging contributions in the thermo- 
dynamic limit is physically unreasonable. The prob- 
ability of a trajectory {X{t),X{t)} is given by P oc 
/ V[X,X]e-^^^^'^\ with action S[X,X] ^ f dt[XX - 
X{aX + bX^ + ...)+ X^{a'X + b'X"^ + ■■■)]■ As long 
as the network has finite average degree (7 > 2), the 
action S should be finite, whereas the diverging coeffi- 
cients contain terms of order Zp ~ N'' ^ (for 
kc ^ N^/^). As we will directly verify later for several 
models, a physical action is recovered when assuming the 
existence of some X-dependent upper cut-off k*{X) such 



that above this value the dependence of Xk on k is weak 
and does not generate divergences. In other words, close 
to the critical points, the system should adjust itself in 
order to eliminate unphysical divergences, with a mecha- 
nism similar to the one assumed by Goltsev et al. in their 
phenomenological theory of equilibrium phase transitions 
on networks [5|. 

In the following, the validity of this general method 
is demonstrated discussing three applications to relevant 
non-equilibrium processes on scale-free networks. 

Voter Model - In the voter model (VM) each 
node of the system is endowed with a binary variable 
(opinion). At each time step, a randomly chosen node 
selects a neighbor at random and assumes its opinion. 
If wc call qk{t) the probability that a node of degree k 
has opinion f , the corresponding rate equation (HMF) is 
qk{t) = Qit) - qk{t), where Q{t) = Y.kkP{k)qk{t) / zi is 
the probability to select a neighboring node with opin- 
ion 1. According to the previously exposed theory, we 
identify X{t) = Q[t) as a good descriptor of the global 
dynamics of the system. Mean and variance of the varia- 
tions AQ can be easily computed on uncorrelated random 
graphs in the MF approximation [l^. Q is a conserved 
quantity, so E[AQ] = and we arrive to the fluctuations- 
driven Langevin equation 

m = ^^Qit) [1 - QWMt) (1) 

where ri{t) is a delta correlated white noise. Processes 
with a global conserved quantity have a trivial infinite- 
size limit, but it is still interesting to study the effects of 
the degree heterogeneity on the finite-size scaling of r, the 
characteristic time needed to reach an absorbing state. 
Here Q(t) performs a kind of random- walk and reaches 
one of the absorbing states {Q = 0, 1) in a characteristic 
time T cx Nzf/z2- For a power-law network with cutoff 
kc iV^/", we have t - TV for 7 > 3, r - TV/ logiV for 
7 = 3, and T ~ n{'^~3+j)/u^ for 2 < 7 < 3, in agree- 
ment with the results by Sood and Redner Q. Studying 
directly the stochastic evolution of q{t) = J^k P(k)qk{t) 
would have led incorrectly to the mean-field result t oc N 
also for 7 < 3. 

Our approach can be applied to the reverse VM (or inva- 
sion process) as well, for which the degree-dependent rate 
equations suggest to study the global quantity R{t) — 
J2k P{k)ziqk{t)/k. As previously Q{t), now R{t) is a con- 
served quantity, and the corresponding Langevin equa- 
tion gives T (X N 

Glauber-Ising dynamics - With respect to the Voter 
model, the Glauber-Ising dynamics has non-zero drift 
and allows to directly verify the relation between our 
method and the equilibrium theory of critical phenom- 
ena on networks. Consider the finite-temperature dy- 
namics on a general graph [18]: on each node i, there 
is a binary spin cr,; = ±1 and spin-flipping events oc- 
cur with rate [I -t- tanh(/3 ^^.^^^-^ ""jOl/S {v[i) are the 
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neighbors of i). In the annealed network approximation, 
one can easily derive the rate equation for the degree- 
dependent magnetization rhk — —m^. + tanh(/?A;/i(t)), 
where h — P{k)hk — '^i^kP{k)mk/ zi is the aver- 
age local field. Note that in the quasi-static approxima- 
tion mfc — tanh(/3A:/i), that close to the critical point can 
be approximated by oc kh for k <^ k* — \/[3h and 
mfc « 1 for /c ^ fc* = l/(3h. Hence, hk oc k'^h/zi for 
k <^k* and hk ~ I for k:s> k*. 

A detailed study of the stochastic evolution equation for 
h{t) will be presented elsewhere Here we consider 

the mean-field rate equation surviving in the infinite size 
limit, i.e. 

^h{t)^-h{t) + Y^ tanh [f3kh{t)] . (2) 

fc 

Expanding the hyperbolic tangent for small h ^ 1, we 
have h ~ -h + Bh^ ~ B'^h?'^. The behavior of the 

z\ ' 3zi 

system depends on the relation between the real cut-off 
/cc and the dynamic cut-off k*(t) (x l/h{t). At the critical 
point h{t) relaxes to zero, therefore there is a typical time 
tc such that k*{tc) — kc, after that any finite system be- 
haves in a mean-field way, except for some anomalous A^- 
dependent prefactors (and the effect of the noise that we 
have neglected [l^ ) . This is true also for infinite systems 
with finite maximum degree {kc < oo). In these systems 
the phase transition occurs at /3 = /3c = Zi/z2 (x k2~^ 
The critical relaxation dynamics m{t) ~ is 
given by solving h oc —^h^ and assuming that close to 

the critical point m = P{k)mk — h. The standard 
MF exponent z = 2 is recovered, but there is a cut-off 
dependent prefactor that accelerates the dynamics. 
On infinite networks with diverging degree cut-off, the 
sums over k have to be performed with the modified up- 
per bound k* oc 1/h to avoid unphysical divergences. In- 
deed for much larger degrees the Taylor expansion itself is 
not correct, because the hyperbolic tangent is close to 1. 
In other words, g{k, h) = tanh(/3fc/i) for the Glauber dy- 
namics. Performing the averages with cut-off k* and ne- 
glecting higher orders in h, we find h{t) ~ —h + A{/3h)"'~'^ 
for 2 < 7 < 3 is a constant), whereas the second 
term gets a logarithmic correction ~ Aphlog{const/h) 
for 7 = 3. There is no paramagnetic phase at finite 
temperature {(Be = 0) [1] and the critical dynamics is ex- 
ponentially fast. For 3 < 7 < 5, the second moment is 
finite but (fc*) — > 00, therefore the effective theory reads 
h(t) = -h + Al3h - B{(3h)^-'^. The critical tempera- 
ture is now finite and the critical dynamics presents the 
anomalous MF exponent z = 7 — 3. The standard MF 
exponent z = 2 is recovered only for 7 > 5. 
Our results corroborate the idea previously suggested 
in Ref.|J] that the non-equilibrium Ising model on net- 
works could be studied using the dissipative dynamics 
dm/dt = -^^^ of soft spins under a generalized 7- 
dependent Landau- Ginzburg potential J-^. 



Reaction- Diffusion processes - In reaction-diffusion 
processes [2^], particles are free to diffuse on a network 
and interact with a given rate when they are on the same 
site. We are interested in "bosonic" process, i.e. there 
are no restrictions on node's occupancy [13] ■ 
The natural degree-dependent quantity to study is the 
average number of particles pk in nodes of degree fc. 
Each variable satisfies the rate equation pkif) = —pk + 
kp/zi + A[pfc], where the first two terms are due to dif- 
fusion in uncorrelated random graphs and the third is a 
reaction kernel that depends only on pk- For simplicity, 
we focus on the dynamics of branching-annihilating par- 
ticles, i.e. 2A and A 2A, with reaction kernel 
A[pfc] = ppk - Ap2. 

The global variable contained in the degree-dependent 
rate equations is the particle density p. Its relation 
with Pk is obtained self-consistently by solving the rate 
equations in the quasi-static approximation, Pk{t) ~ 
i^[l - ^1 + 4Afcp(i)/((/i- l)2zi)]. Now that we have 
Pk = g{k,p{t)), we develop g for small values of p, close 
to the critical point fic = 0. We find two regimes: pk{t) oc 
kp{t)/zi for fc ^ fc* cx zi/p{t) and pk{t) oc ^Jkp{t)/ zi for 
fc > fc*. 

Starting from the degree-dependent variations Ap^ and 
averaging over P{k), we find the Langevin equation 

p{t) ^^ip- X{pl) + A-i/y/ip + 2X{pl)fj{t). (3) 

According to the above relation between pk and p, the 
average over P(fc) can be splitted in two integrals: one 
for k < k* and the other for k > k*. In a finite network 
with maximum degree fcc, 

f 2 77 { p^^^ + const X k'^~^ p, kc > k* 

J P{k)pldk0C^^ fc3-V, kc<k*. 

Close to the critical point, the particle density is small, 
thus we expect that on finite networks and infinite net- 
works with finite cut-off the system will eventually en- 
ter the mean-field regime (fc* > fcc). The corresponding 
Langevin equation is 

p{t) ~ap^ h'p" + ^{ap + b'p^)/N7]{t) (5) 

in which the various constants have been absorbed in two 
effective reaction rates a oc p, and b' = 2hk^~^ {h oc X). 
Assuming fcc ~ A^^/", the stationary density in the active 
phase follows a mean-field scaling pact oc A{N)p, with a 

non-trivial prefactor A{N) A^ ^ . For a generalized 

1 _ 9+1-7 

reaction qA 0, we have pact oc pi-^ N . For fi- 

nite A^ the active state is metastable, that is the system 
will eventually decay to the absorbing state because of 
noise- induced fiuctuations. The amplitude of the lead- 
ing noise term scales like A^ ^ . In case of networks 
with structural cut-off (w = 2), the noise is always small 
at large but finite A^. If we prepare a system in the 
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metastable state {pact 1/vN) just above the critical 
point (p > 0), it will decay to the absorbing state in a 
time e^(^''' \ But in a network with natural cut-ofF 

= 7 — 1 , the dynamical fluctuations become patholog- 
ical in the limit 7-^2: the density of the metastable 
state is just 0{1/N) and it will decay to p = in a finite 
timescale. The problem of large fluctuations induced by 
multiplicative noise in scale-free networks requires fur- 
thcrly investigations [3]. 

In infinite networks with diverging degree cut-ofF kc, at 
every time we have k*{t) < kc, so using Eq|4]we get the 
effective rate equation p(t) ~ ap — bp'^~'^. The exponent 
7 — 1 is universal, independent of the number q of par- 
ticles involved in the annihilation reaction. In summary, 
the general effective mean-field equation for the critical 
behavior of reaction-diffusion processes on uncorrelated 
random networks with P{k) ~ k~'^ (7 > 2) is 



p{t) ^ ap — bp 



(6) 



where the second term dominates over the third for 
7 — 1 < 5 and viceversa. Logarithmic corrections appear 
for 7 — 1 = 9- This expression for two-particles anni- 
hilation, q = 2, resembles the phenomenological theory 
proposed by Hong et al. |14| . However, our expression 



is strictly correct only in the infinite-size limit, when the 
stochastic term due to the granularity of the system is 
absent. The addition of a multiplicative noise term, such 
as ^/prj proposed in [l3 |. wherever justified, should not 
change the mean-field behavior obtained from Eq. [6l 

A similar analysis could be performed for "fermionic" 
(or hard-core) particles as well (i.e. nodes are occupied by 
at most one particle) [1, [l3|- For instance, the fermionic 
branching-annihilation model is described by a rate equa- 
tion of the form pk{t) = —pk{t) + kp/zi + A[pk,p], with 
HPk,p\ = i^kil - Pk{t)]p{t)/zi - 2kp{t)pk{t)/zi {v is an 
effective reaction rate) [10]. At fixed global density, we 
observe the scaling pk{t) oc kp{t)/zi for fc <C fc* cx zi/p{t) 
and pk{t) ~ const. < 1 otherwise. Hence previous calcu- 
lations can be easily extended to fermionic models, ob- 
taining qualitatively similar results [l9| . In particular, 
our analysis corroborates recent results by Boguha at al. 
[1^ for the contact-process on networks. 

Conclusions - We have proposed a general method to 
derive improved mean-field theories for non-equilibrium 
critical phenomena on heterogeneous networks. It al- 
lows to choose the correct mean-field order parameter for 
many different models, and obtain the anomalous mean- 
field scaling laws in a systematic way for both finite- 
size and infinite-size systems. This approach is simpler 
than the widely used HMF method and represents the 
non-equilibrium generalization of the equilibrium theory 
of critical phenomena on scale-free networks [5]. The 
method was successfully applied to the voter model, the 



Glauber dynamics and reaction-diffusion processes. For 
all of them, the nature of critical behavior depends on 
the cut-off of the degree distribution. For finite cut-offs, 
the system eventually recovers mean-field behavior close 
to the critical point, whereas fluctuations may strongly 
affect single realizations of the process for very heteroge- 
neous networks [l9j . On the contrary, infinite networks 
with diverging degree cut-off present an anomalous mean- 
field critical behavior, directly dependent on the expo- 
nent 7 of the degree distribution. 
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